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ABSTRACT 



Experiments were conducted at the NURDC oceanographic tower off San 
Diego, California to describe the small scale physical properties in the 
upper ocean and to determine their temporal and spatial interrelation- 
ships. The measured parameters included water particle velocities, tem- 
peratures, salinity, sound speed, sound speed (phase) modulation, sound 
amplitude modulation and surface waves. Simultaneous time series meas- 
urements are analyzed for a twenty minute record taken at a depth of 
7.3 m. The spatial correlation distance for temperature was found to be 
about 25 cm, Sound modulation of a 60 kHz sound beam underwent a mean 
attenuation of about 57o. The correlation distance of an effective in- 
homogeniety due to the combined temperature, salinity, and bubble effects 
was found to be 8.8 cm. The change in the index of refraction of the 
speed of sound at 60 kHz (a bubble resonant region) was found to be 99% 
induced by bubbles. The energy density spectra were dominated by wave 
induced motion. First order wave theory accurately predicts orbital 
water particle motion and its irrotationality . High coherence is found 
between wave-induced water particle motion and sound modulation and de- 
creasing coherence between water particle motion and salinity, sound vel- 
ocity, and temperature. 
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INTRODUCTION 



A. BACKGROUND 

The importance of detailed studies of sound propagation in the upper 
layer of the ocean has become increasingly apparent as the role of the 
submarine continues to grow in strategic and offensive capabilities. Also 
significant from a scientific viewpoint is that the manner in which sound 
propagates can be used to effectively describe properties of the medium. 
Thus in order to understand sound propagation, a detailed knowledge of 
the parameters that affect it and their interrelationships with each other 
and sound is required. Wilson [ lO] developed an equation that related 
the speed of sound in water to temperature, salinity, and pressure. Lie- 
bermann [4] delved into the effects of temperature inhomogeneities on 
sound propagation and determined that inhomogeneities do have an effect 
on the intensity and refraction of propagated sound waves. Mintzer and 
Stone [&3 theoretically developed an equation that related to the reduction 
in sound intensity of a propagated wave to wave number, correlation dis- 
tance , distance traveled and changes in the refractive index of the medium. 

In November of 1971, Bordy [2], Duchock [s], Rautmann [?], Seymour 
LS], and Smith [9] from the Naval Postgraduate School took measurements 
at the NURDC tower of various oceanographic parameters in the near surface 
zone while propagating sound waves over prescribed paths. The parameters 
measured included sound velocity, sound amplitude and phase modulation, 
water particle velocities, temperature, salinity, and wave height. The 
main purpose of this paper is to extend the results of the references just 
cited carrying out a more detailed study of these paramenters and their 
interrelationships . 
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B, OBJECTIVES 



1. To determine the interrelationships between the various para- 
meters (temperature, salinity, sound velocity, modulation and wave height) 
in the near surface ocean zone. Only those processes v^;ith periodicity 
between 0,5 to 50 seconds will be considered. 

2. To compare results with established equations and relationships. 
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II. EXPERIMENT 



On 21 October 1971, Bordy, Duchock, Rautmann, Seymour, and Smith 
from the Naval Postgraduate School conducted experiments at the Naval 
Undersea Research and Development Center Oceanographic Research Tower 
at San Diego, California. The purpose of the experiment was to describe 
the small scale physical properties in the upper ocean and to determine 
their temporal and spatial interrelationships. The tower is approximately 
one mile off Mission Beach in a mean water depth of 17 meters. The 
structure is made of steel and concrete and built into the sea floor. An 
enclosed upper deck provides space for electronic equipment while an open 
lower deck provides a platform to place instruments in the water. On the 
northern, southern, and western sides are located winch controlled in- 
strument carriages that can be adjusted to any desired depth. (Figure 1 
and Figure 2) . 

The parameters measured during the experiment included orthogonal 
water particle velocities, an array of temperatures, salinity, sound 
speed, sound speed (phase) modulations, sound amplitude modulation and 
surface waves measured directly over the column of interest. All sensors 
were placed within a 2 meter cube of water. Simultaneous time series 
measurements were made for twenty-minute durations. The experiments 
were conducted at four depths. In this paper. Run Five, at a depth of 
7.3 meters, will be analyzed in detail. The spatial configuration of 
the sensors is shown in Figure 3. The actual geometrical separations 
are shown in Figure 4. 

During the 24 hour period of the experiment, a bi-modal surface 
swell of about 0.3 meters height was visually observed with one component 
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Navy Undersea Research and Development Center 
Oceanographic Research Tower 

Figure 1. 
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Research Tower Schematic 
Figure 2. 
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Figure 3. Spatial Configuration of Instruments (from Seymour [8]) 
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having a period of about 7 seconds from 285^ true and the other compo- 
nent with a period of about 5 seconds from 260^ true. The surface was 
observed to be non-breaking with a wind capillary system present. The 
wind was from 310^ true at a speed of about 8 knots. A moderate temper- 
ature gradient was present with a temperature of 17.5^ C at the surface 
decreasing to 15^ C at a depth of 18 meters. 
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III. INSTRUMENTATION 



A brief description of the instruments used will be given in this 
section. A more detailed description can be found in the papers of 
Bordy, Smith, Duchock, Rautmann, and Seymour. Instrument details are 
given in manufacturers* technical manuals. 

A. SOUND VELOCITY SENSOR 

Sound velocity was measured with a Ramsay MK-1 deep sea probe util- 
izing only the sound velociraeter portion of the probe. The Ramsay veloc 
imeter is of the "sing-around" type consisting of a transducer and 
transmitter which transmits a pulse of 4 MHz into a 25 centimeter sound 
path. The pulse is reflected twice to reduce errors due to water motion 

B. TEMPERATURE SENSOR 

Temperature was measured using a specially constructed Wheatstone 
bridge circuit utilizing a temperature ’ sensitive thermistor as one leg 
of the circuit. For Run Five, two such circuits were employed. The 
type of temperature sensor used was a type K 496 isocurve thermistor 
manufactured by Fenwall Electronics, Inc., Farmington, Massachusetts. 

The specific thermistors were extracted from an expendable bathythermo- 
graph. It is a double-bead thermistor which incorporates a special 
aging process, giving stable and repeatable accuracies not usually ob- 
tained with ordinary thermistors. These thermistors have a time con- 
stant of 150 milliseconds. 

C. BAYLOR WAVE GAUGE 

The Baylor Model 13528 is a wave profile recording system which 
measures instantaneous sea surface elevation as a function of time. 
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Measurements are accurate within one percent of actual height. The sys- 
tem consists of three components; the wave staff, transducer and recorder. 
The transducer unit measures the electrical resistance of the length of 
the staff above the water surface. The conducting seawater provides a 
short circuit across the cables (wave staff) immersed in the water. A 
direct current voltage proportional to the height of the sea surface is 
taken across the resulting fluctuating resistance. 

D. WATER CURRENT METER 

The water particle velocities were measured with an Engineering 
Physics Company water current meter Model EMCM-3B. The EMCM-3B is an 
electromagnetic current measuring device that measures two orthogonal 
components of water particle velocity from zero to five meters per sec- 
ond. Two d.c. voltage signals proportional to the magnitude of the 
orthogonal velocity components are produced. The measured response is 
flat to at least 0.5 seconds, the high frequency limit of the calibra- 
tion facility. The maximum output error is one percent of full scale 
reading. 

E. SALINITY SENSOR 

The salinity sensor employed was part of the Bisse t t-Berman STD 
Model 9006. An inductively-coupled sensor detects conductivity, which 
is also compensated for temperature and pressure effects, thereby pro- 
ducing an output totally dependent upon salinity. The time constant for 
the salinity sensor is .35 seconds. 

F. MODULATION 

In order to measure environmentally caused amplitude modulation in 
the ocean, a high frequency acoustic carrier wave was directed at a 
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receiving hydrophone mounted along the acoustic axis at a separation of 
two meters. An envelope detection system was used to measure the ampli- 
tude modulation of the received signal, A USRD type F-27 transducer and 
an Atlantic Research Corporation type LC-10 hydrophone were used in the 
experiments , 

Go SOUND VELOCITY (PHASE) 

In determining sound velocity by phase measurements, one F-27 trans- 
ducer was used along with one LC-10 and one LC-5 hydrophone. See Rautmann 
[ 7 ] for further details. 
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IV. DATA PRE-PROCESSING 



Run Five was started at 1728 and ended at 1748 on 21 October 1972. 
During the run, data was recorded on a Sangamo Model 3500, 14 channel 
FM tape recorder. The data consisted of temperatures 1 and 2, salinity, 
wave height, u and w water particle velocity components, sound modula- 
tion, and sound velocity. The taped data was transcribed to a recti- 
linear eight channel Clevite Brush Recorder after amplification by 
Krohn-Hite Model 3322 variable filters. 

A Calma Company Model 480 mechanical digitizer was used to digitize 
the continuous strip chart information at 100 points to the inch. This 
corresponds to a sampling interval of 0.2032 seconds. This information 
was recorded onto a 7 track magnetic tape which is compatible with the 
CDC 6500 computer at the Fleet Numerical Weather Central, Monterey. The 
data was transferred to cards using program C0 N\TlRT (Lynch, 1970) so that 
it could be read for analysis by the IBM 360 at the Naval Postgraduate 
School which uses 9 track magnetic tapes. The calibration factors used 
to convert the data into the appropriate units are listed in Table I, 
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TABLE I 



CALIBRATION FACTORS FOR DATA 



CALIBRATION FACTOR 


MEASUREMENT 


UNITS 


.23 


Temperature 




.508 


Salinity 


ppt 


00 


u , w 


meters/second 


.34 


Sound Velocity 


meters/second 


.211 


Fractional Modulation 




.775 


Wave Height 


meters 


This information is 


included here to facilitate 


future analysis of 



data . 
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V. ANALYSIS OF DATA 



A. FOURIER FILTER 

The experiment was designed to study processes in the approximate 
frequency band 0.02 to 2 Hz corresponding to time scales of 0.5 seconds 
to 50 seconds. The initial analysis of the raw data showed that the 
energy density spectra of some parameters had excessive leakage from 
the very low frequencies. The time correlations were of the order of 
minutes. Any effects caused by forcing functions outside the frequency 
band of interest, such as tidal or internal wave effects, will not be 
considered. Thus it is necessary to eliminate the low frequency com- 
ponents in the raw data records. This was accomplished by using a 
Fourier Filter. 

The first step is to generate a Fourier series of the record. The 
fast Fourier Transform algorithm is used to conserve computer time: 

CO 

F(l) = ^ y (A cos net + B sin njut) 
z n n 

n=l 

where F = parameter 

t = time 

n'jj = nth harmonic of fundamental radial frequency it 

Next, a time series is generated using only the low frequency components 
that are to be eliminated. This time series, composed only of the low 
frequency components, is then subtracted from the original time series 
for each value, and the remainder is the filtered value of the parameter 
at each given time. In this experiment, the first fifteen low frequency 
components were used to generate the filtered time record, corresponding 
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to the components with a fifty-five second periodicity and greater. Al- 
though certain records needed filtering vzhile others did not, all records 
were filtered for consistency. 

The effect of the Fourier Filter on the temperature 2 spectra can be 
seen in Figure 5. Notice the removal of the energy from 0 to .018 Hz in 
the energy density spectrum. After this point, the energy spectra coin- 
cide for the remainder of the spectrum. Also of importance is that the 
phase and coherence spectra remain unaltered in the region of interest, 
greater than .018 Hz. 

B. SPECTRAL ANALYSIS 

The energy-density spectrum is calculated from the auto-covariance 
function of the recorded time series: 



where the auto-covariance function is expressed as a function of the lag 
time, T. If the auto-covariance is normalized by division by the vari- 
ance, it is called the auto-correlation function. The Fourier transform 
is the power’ or energy-density spectrum. 



A Parzen lag window was applied to the functions to account for the 
finite length of the records in computing the spectra. Since the Parzen 
window has no negative side lobes, this eliminates any numerical insta- 
bility problems that could arise in the analysis. 

Let P(t) denote the Parzen Window Function, then the modified energy 
density spectrum can be written: 




dt 




-i2nfr 



dr 
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5ll (f) = 


= J P(t)0ii (T)e dT 

- 00 



The energy-density spectrum is the statistical estimate of the energy 
density distributed over the frequency bands of the spectrum. The total 
potential energy-density of the sea surface height is proportional to the 
total area under the wave height spectrum. Similarly, the area under 



the velocity spectrum is 


proportional to the kinetic energy-density. The 


area under each curve is 


equal to twice the variance of the random vari- 


able . 





The cross-spectral density is defined as the Fourier transform of 
the cross-covariance function. The cross-covariance function is similar 
to the auto-covariance function except that a second time series record 
is lagged against the first instead of a single record lagged against 
itself. The cross-covariance function indicates the dependence the 
values of one record have on that of another for various lag times. 

The cross-covariance function is: 



01 2 (T ) = 


1 

= lim - J X (t)x (t + t) dt 

T--® -T/2 



The modified cross spectral density is 



§12 (f) = 


= J P(t)0i2(t) e dT 

.00 



The cross-spectrum can be defined in terms of its real and imaginary 
parts : 

^12 (f) = Ci2 (f ) “ iQisCf) 

where the real part of ^12 (f) is known as the co-spectrum: 



Cl2(f) = 


00 

= 2 J* L^i 2 (t) + §12 (-t)] cos (2nfr) dr 
0 
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■T 







a 



and the imaginary part of ^12 (f) is the quadrature spectrum: 

00 

Qi2 (f) = 2J [§12 (t) 4- 9i2(~t) 3 sin (2TTfT) dT 
0 

The spectral phase angle represents the average angular difference 
by which the cross correlated components of lead those of x^(t) 

in each spectral frequency band. The phase angle is defined by: 



fi'is (f) = arctan 



Lci2(f)J 



(f) ranges from - 18 CP to -fl 80 P . 

The coherence function is defined by: 



Yii(f) 



[^ir (f)]^ 

•11 (^)^2 2 (^) 



with the provision that: 

I ^12(f)| ^ ^11 (f)^2s(f) 



which results in: 

0 ^ yJ (f) ^ 1 



The coherency spectrum is a measure of the correlation between two time 
series as a function of frequency. Three possible causes of coherency 
values being less than unity are: th^ system is not linear, excess noise 

is present in the measurement, or the output time series is due to the 
input time series as well as other functions (Bendat and Piersol Ll]). 



C . SPATIAL CORRELATION 

The data were measured as a time series record at specific points in 
space. The majority of theory in this area of study requires knowledge 
of spatial correlations. Hence, meaningful comparisons cannot be made 
without measures of spatial correlations. The operation described below 
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is designed to convert the time domain record to a spatial one. The tem- 
poral measurements are transformed to a Lagrangian frame knowing the ve- 
locity and acceleration at the point of parameter measurement. 

It is assumed a priori that the correlation lengths are small com- 
pared to the spatial scales of motion. At 7.3 meters depth the orbital 
excursion of a 10 second wave having a rms wave height of 30 centimeters 
(mean conditions during experiment) is 80 centimeters. It is further 
assumed that the temperature field is isotropic to simplify the analysis. 
Measurements of velocity in the X and Z plane were taken simultaneously 
along with temperature. Although the velocity components were measured 
a small distance from the thermistor it is consistent to assume that the 
velocity components can be transferred to the thermistor location, r^, 
due to the small separation distance involved. If one now considers the 
thermistor moving with the velocity, U, and acceleration components, it 
will travel a path in the X-Z plane. At each instant the velocity is 
measured, the distance traveled by the thermistor, r, is calculated by 



, AU^t . 

r = r -{ + UAt 

o 2At 

(At= time increment) 



(C-l) 



whereupon the next temperature is assigned to this distance ”r.*' Thus 
an array is created where a temperature is assigned for a certain dis- 
tance *'r” in the X-Z plane. For each temperature measurement in time 
there is an associated distance giving T(r,t). In order to facilitate 
calculation of the auto-correlation function, any holes in the array 
were filled by linearly averaging the temperature over the distance 
increment . 

The question arises as to the validity of the concept in relation to 
the occurence of a negative velocity which would necessitate a backward 
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movement of the thermistor, bringing about the possibility of two differ- 
ent temperatures being assigned to the same point in the X-Z plane. The 
probability of this happening is small, and the effect of its occurrence 
is negated by averaging the first temperature with the second one. If 
by chance a third temperature were assigned this point, it would be aver- 
aged with the previously averaged temperature; hence the last temperature 
would alv;ays have a weighting factor of 50 per cent. 

A similar scheme using an Eulerian frame can be derived by expanding 
the temporal measurement at a point in a Taylor series: 



T(r t) = T(r t) + + 



(C-2) 



v;here Ar = r - r^ and Ar is small. Employing Taylor's hypothesis that 

Ar 

over short distances and time U = ~ , Equation (C-2) can be rewritten: 



T(r t) = T(r t) 4- + 

Ur,c; i'.rQ.c; + 2U^At^ 



(C-3) 



The two schemes are equivalent, but computationally the Lagrangian form 
is more convenient and is employed here. 



D. AUTO-CORRELATION FUNCTION 

With the conversion of a time series record to a spatial record in 
the preceding section, the description of the parameter in the medium by 
the spatial auto-correlation function is now possible. The importance of 
the spatial distribution of the parameter is illustrated in the following 
example. A sharp temperature boundary occuring within an acoustical 
wavelength of the propagated sound may cause reflection. The linear 
dimensions of the temperature inhomogegeity may cause scatter, reflection, 
or focusing of sound waves. The physical features of a parameter may be 
characterized by the auto-correlation function. 
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The compulation of the spatial auto-correlation function, 0(p,t), 



from the temperature record, T(r,t) is defined: 



. L/2 

(j5(p,t)=-J T(r,t)'T(r+p,L) dr (D-1) 

^ -L/2 

0(0, t) 



where p is the lag distance and L, the tota] length of the field. Fig- 
ure 6 shows the spatial auto-correlation function for temperature 1 and 
2. Figure 7 sho\^s the time auto-correlation function for the same para- 
meters. Note the rapid decay to zero by both functions and also the 
sinusoidal components; thus both cases reveal the presence of periodic 
components, namely wave induced particle motions. 

The distribution of temperature 2 in time is shown to closely ap- 
proximate a Gaussian distribution in Figure 8, Figure 9 (using the 
Fourier Filter). Thus a good description for the function over short 
time lags would be: 

0(t) = e" (D-2) 

where t^ is the correlation time. Assuming a similar spatial Gaussian 
distribution gives: 



6(p) = e 




(D-3) 



The mean **size" of the inhomogeneities might reasonably be described 
by the value of p at which the temperature 0(p) drops to 1/e denoted by 
po in Equation D-3. These distances are known as correlation distances 
(Lieberraann [4]). For temperature records 1 and 2, their values of Po 
were 26 and 27 cm respectively. 

The consistent results for correlation distances show that the con- 
version from time to the spatial domain is reasonable. Hence, the 
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FIGURE 6 . SPATIAL TEMPERATURE AUTOCORRELATION FUNCTIONS 










i 






1.0 




TEKAPERATURE 1 




FIGURE? .TIME TEMPERATURE AUTOCORRELATION FUNCTIONS 



31 




#1 




1 




FIGURES. TEMPERATURE 2 PROBABILITY DENSITY FUNCTION 
WITHOUT FOURIER FILTER 




Figure 9. Temperature 2 Probability Density Function With Fourier Filter. 
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approximate mean size of the temperature inhomogeneities at a depth of 
7.3 meters are characterized at about 26 centimeters. The difference in 
correlation times for temperature 1 and 2 as seen in Figure 7 (4.5 sec- 
onds versus 1.4 seconds) could be due to higher turbulence in the loca- 
tion of thermistor 2 which would result in smaller inhomogeneities. An 
alternative possibility would be a malfunction of the thermistor. Regard- 
less, the conversion of a time series to a spatial domain has the effect 
of averaging out the localized turbulence by filling the spatial tempera- 
ture array with smoothed values of temperature variations. Thus the 
spatial processing tends to average by minimizing the effect of localized 
turbulence . 

A satisfactory distribution function for the near surface zone must 
take into account periodic motion induced by surface waves. Therefore, 
a more accurate description of the auto-correlation function would be: 

0(P>t) = j~ e *^n-l^^n^°^ n(kx-ot) + n(kx-cyt)]J 

- ^ 

However, this function is entirely too cumbersome. R(p) = e Pq will be 
used because it describes the function very well through the correlation 
distance up to the first zero crossing and is capable of being integrated 
which is important in a later section. 

E . WILSON ' S EQUATION 

Wilson *s Equation [2] was originally developed to relate the speed 
of sound to the properties in ocean water that affect it, namely pressure, 
temperature, and salinity. An abbreviated form of Wilson* s Equation is: 

c = 1449 + 4.62T - 0.055T^ + 0.00031^ + (1. 39- . 012T) (s-35) + 0.017z 
where c “ speed of sound in meters per second 

T ~ temperature in degrees Centigrade 
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s - salinity in parts per thousand 
z = depth in meters 

One of the initial aims of the experiment was to show that the vari- 
ance of the speed of sound was equal to the sum of the variances of its 
constituents multiplied by their respective coefficients. This would 
establish the validity of Wilson* s Equation in terms of variances as 
well as indicate the accuracy of the temperature, salinity, and depth 
data. 

Wilson* s Equation applied in this manner was not successful for the 
experiment. Let 

var c = ( (c-cf) 
dT = < T-T > 
ds = ( s-s ) 

where the overbar implies averaging and ( ) indicates an ensemble average. 
After simplifying by retaining only first order terms of temperature and 
salinity and dropping the small depth term, 

var c = ((5.04dT + 1.39 ds)^) 
var c = ((25.40dT^ -f 14.00 dT*ds + 1.93ds^)) 
and dT'ds = covariance (T*s) = cov (T*s). 

In order to measure the variance of c using Wilson* s Equation, both 
temperature and salinity must be measured at the same point at the same 
time. The time constant for the salinity sensor was .35 seconds, and 
for the thermistors .15 seconds. The average spatial distance between 
the instruments was 2.14 meters. 

Correlation distance has been described in a previous section as 
representative of the mean size of the inhomogeneiety . In order to en- 
sure that temperature and salinity sensors provide accurate results. 
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they must be located v/ithin the same inhomogenci ty along with the sound 
velocity meter. This is obtained only if all the instruments arc placed 
within the smallest correlation distance of the parameters involved. As- 
suming that the mean size of a salinity inhomogeneity would be no smaller 
than that for temperature, then all instruments v/ould have to have been 
placed within 26 centimeters in this experiment (as previously determined 
from temperature correlation distances) in order to provide an accurate 
evaluation of Wilson's Equation. For the experimental data taken, an 
application of Wilson's Equation in this case would be futile as Seymour 
[S] showed in his computer results. 

As the distance between instruments increases in the random medium, 
one would expect that local mechanisms operating within one correlation 
distance to bear less and less resemblance to mechanisms operating at 
increased correlation distances; hence, the temperature and salinity 
would approach statistical independence with increasing distance, causing 
cov (T*s) to approach zero by definition. This was found to be the 
case in Run Five where the cov (T*s) integrated over the spectrum was 

found to be -.8 x 10 ^(C^-PPT). The variance of temperature was 1.6 x 

o 2 2 2 

10 (C ) , the variance of salinity 0.22 x 10 (PPT) , and the variance 

“A- 2 

of sound velocity 1.1 x 10 (raeters/second) . 

The importance of the waves to the variance of the parameters is 
indicated by band-pass filtering. Band-pass filtering was accomplished 
by high-pass filtering with the Fourier filter described before and 
calculating the variance by summing the spectral contributions only of 
frequencies up to 0.16 Hz. The band-pass filter was chosen because it 
included that part of the variance where organized wave-induced particle 
motion was dominant in its contribution to the variance, and the effects 
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of turbulence and instrument noise were insignificant. Table II shows the 



results of the band-pass filter. 



TABLE II 

VARIANCES OF SOUND VELOCITY, TEMPERATURE, AND SALINITY 

Variance Entire Spectrum Band-pass Filter % Variance 

in wave dom- 
inated region 



var 


c 


1.1 


X 


lO"'^ 


1.0 


X 


lO"^ 


91% 


var 


T 


1.6 


X 


lo"'^ 


0.6 


X 


o 

1 


38% 


var 


s 


21.5 


X 


o 

1 


19.1 


X 


lO"'^ 


90% 



On inspection, the high value of variance for salinity makes it 
questionable as one would expect a value of the same order of magnitude 
as the temperature, not 13 times as great. Further investigation in- 
dicates that the measured conductivity is compensated continuously for 
pressure and temperature changes with the result that the indicated sal- 
inity has a tendency to oscillate; thus the high value for the variance 
of salinity should be disregarded. This is significant in itself, and 
when combined with a lack of the covariance term, cov (T*s), it makes 
any proof of Wilson’s Equation only conjecture. Hence, future investi- 
gations must use more appropriate salinity instrumentation in the near 
surface region of the ocean. 

The frequency of the sound velocimeter was 4 MHz. In this frequency 
range, the bubbles present in the near surface zone have no effect on sound 
propagation due to the extremely small wavelength of the sound wave. 

Thus by measuring the variance of the sound velocimeter and comparing it 
by Wilson’s Equation to the variances of the temperature , salinity and 
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their respective covariance, one gets a crosscheck on the refraction of 
sound in the medium due only to salinity and temperature. 



F. AMPLITUDE MODULATION IN A PJLNDOM MEDIUM 



In an ideal homogeneous undisturbed medium, a sound wave propagated 
at a fixed frequency and amplitude arrives in the far field undiminished 
except for spherical spreading and absorption losses. An adequate de- 
scription of the process for a homogeneous medium is given by the general 
wave equation: 



, , 2 d^P 
V P = 1/c — 5 



(F-1) 



where V = LaPlacian operator 

P = pressure 

c == speed of propagation in medium. 

The medium under study is the near surface ocean. Bubbles, tempera- 
ture and salinity patches, and scatterers are all present in a random 
distribution. Compounding the problem is the fact that the inhomogene- 
ities are in constant motion due to wave particle motion or convective 
forces. As the sound wave moves through a medium, the inhomogeneities 
cause a variation in the speed of sound which results in refraction. If 
the inhomogeneities were motionless, reflection and refraction of a 
propagating sound beam would cause the amplitude of the sound beam to 
be constant at a given receiver. If the inhomogeneities have freedom 
to move in and out of the path of the sound beam, then the amplitude of 
the sound beam at the receiver will fluctuate. Hence, the wave equa- 
tion must be modified to include spatial and temporal changes in the 
speed of sound. The change in the index of refraction with time and 
space is defined: 
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1 O' 



(F-2) 



„ = c(r,t) ^ Co -f- dc(r,t) ^ ^ clc(r^^ ^ 

Co Co Co 



where: 



Q, ^ .dc,Cr^jJ 
Co 



|i ~ Index of refraction 
co= Average speed of sound 
c = co+dc= Local speed of sound 

Q! = RliS refractive index variation from unity 
t = Time 

r “ Position vector 



Previous measurements of O' for temperature inhomogeneities have 
yielded values such as 1.18 x 10 ^ (Stone and Mintzer in Laboratory [lOj) 
and 7.0 X 10 ^ (Liebermann [b] in Ocean Experiments at depths of 50 
meters). In Run Five of this experiment, the calculated O' for tempera- 
ture inhomogeneities at a depth of 7.3 meters was 3.0 x 10 This was 

calculated by taking the first derivative of Wilson's Equation with re- 
✓ 

spect to temperature and substituting the rms value of the temperature 
variance into only temperature terms up to the third order; the average 
speed of sound was calculated using Wilson's Equation with time averaged 
values for temperature, salinity, and pressure. In comparing Lieber- 
mann *s results to Run Five, the implication is that the temperature in- 
homogeneities are smaller near the surface. This seems reasonable as 
water particle motion caused by wave action decays exponentially with 
depth; hence better mixing would be expected nearer the surface reducing 
the size of the temperature inhomogeneities. 

In their laboratory studies of amplitude modulation in a random med- 
ium, Mintzer and Stone, [bj, defined the coefficient of variation, cv: 

C.v>2 = (F-3) 
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where, cv = The coefficient of variation 

P = The amplitude of the received signal 

The coefficient of variation as defined by Hintzer and Stone is related 
to the fractional modulation measured in Run Five by the following equa- 
tion : 

2 

(.0253) (cv) = var Modulation 

The development of this relationship is found in Appendix I of Smith l9]. 
The majority of theory dealing with propagation in a moderately in- 
homogeneous medium employs the single scattering approximation to the 
wave equation. Mintzer and Stone Lb] have sho\>m that: 

oo 

(cv)^ = a\ f 0(p) dp (F-4) 

o 

where: 

K == Wave number 
o 

b = range from source to receiver 

0(p) = the correlation function taken along the path 

from source to receiver. 

or = root mean square index of refraction variation from 
unity. 

Previously, it was determined that the temperature auto-correlation 

- 

was described by 0(p) = e Po*" • Other parameters affecting the index 
of refraction showed similar Gaussian distributions; therefore, sub- 
stituting the expression for 0(p) (representing a combined temperature, 
salinity and bubble inhomogeneity auto correlation function) into equatiion 
(F-4) and integrating, one gets the following result: 

(cv)^ = ^ /7? K r (F-5) 

o o 

In Run Five the sound modulation was obtained for a sound wave of 
frequency 60 kHz. This frequency has a wavelength of 2.5 centimeters 
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which happens to correspond to a peak in the resonant bubble regime of 

the near surface zone at a depth of 7.3 meters in this experiment, as 

determined from Rautmann L?]. By using two hydrophones to determine 

phase change over a 78 cm path, Rautmann converted the phase change to 

a change in the speed of sound. At 7.3 meters depth during the same 

18 hour time period, an interpolated value for 60 kHz from Rautmann *s 

tables indicates a c(measured) - c (theoretical) = -1 m/sec. C(measured) 

is the speed of sound based on phase measurements while c (theoretical) is 

taken from Wilson *s Equation using average temperature, salinity, and 

depth. Thus the variance of the speed of sound at 60 kHz was determined 
2 

to be 1 (m/sec) due to temperature, salinity, and bubble inhomogeneities. 

Converting this to an a value gives the relatively high value of 
-3 

.67 X 10 . This value is based on bubbles, temperature, and salinity, 

not just temperature alone which gave a value of .25 x 10 Since the 
salinity values were highly suspect and along with the lack of a cov(T’s) 
term, one cannot accurately determine the exact effect of bubbles on O'. 
However if one makes the modest assumption that the salinity contribution to 
dT-ds will be equal to the temperature contribution or less, then the 
temperature and salinity contribute less than one percent of the value 
of a. This shows the enormous importance of bubbles in the near surface 
region upon the refraction of sound at 60 kHz. 

Further proof of the importance of bubbles can be seen in the spec- 
tral analysis of wave height versus temperature, salinity, and sound mod- 
ulation (see Figures 10,11, and 12), In the coherence and phase spectra 
of wave height versus temperature and salinity, small consistent values 
of coherence and phase are found in the wave dominant portion of the 
spectrum (.05 to .22 Hz). However, for wave height versus sound modulation, 
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FIGURE 10 , WAVE HEIGHT AND TEMPERATURE - 2 SPECTRA 



•42 



• ^ 




COHERENCE PHASE SPECTRAL DENSITY (AMPLITUDE -SEC) 





FIGURE 11. V/AVE HEIGHT AND SALINITY SPECTRA 
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FIGURE 12. V/AVE HEIGHT AND SOUND MODULATION SPECTRA 



44 



the coherence value approaches unity and the phase is constant in the 
wave dominant region. This implies that wave height is linearly re- 
lated to sound modulation. Since the effect of the change in pressure 
due to v;ave height is insignificant on the sound modulation at this depth 
then the fluctuation in modulation must be due to the variations in the 
inhomogeneities. With the coherence of wave height and temperature and 
wave height and salinity being low, the process of elimination leaves 
only one type of inhomogeneity -- bubbles. This supports the contention 
for this experiment that the dominating cause of change in the speed of 
sound was due to bubbles as evidenced by very high coherence of wave 
height and modulation. 

This high coherence can be due to two wave related mechanisms: 1) the 

bubbles (inhomogeneities) are moved in and out of the sound path by the 
wave induced motion and, 2) the bubbles are ’’squeezed” by the wave induced 
pressure changes resulting in a change in their resonance characteristics. 

In applying the data from Run Five to Mintzer’s equation, ’’Po/* the 
mean radius of the scattering inhomogeneity, is found to be equal to 8.8 
centimeters. This is not to say that an inhomogeneity exists with a 
radius of 8.8 centimeters, but the effect of bubbles, temperature, and 
salinity provides an ’’effective” radius of 8.8 centimeters for the in- 
homogeneity. 

G. THE EFFECTS OF SURFACE WAVES ON THE MEASURED PARAMETERS 

There are three contributors to the energy density spectrum of hori- 
zontal or vertical particle motion: wave induced particle motion, mean 

convective motion, and turbulence. Mean convective motion was essentially 
zero for the particle velocity in the vertical plane. This leaves wave 
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induced particle motion and turbulence as the two dominant sources in 
supplying energy to the power spectrum of vertical particle motion. 

Turbulence is random three-dimensional motion which can dissipate 
mechanical energy to internal energy through a cascade of eddies of 
diminishing size. The process of turbulent cascade, resulting in dis- 
sipation, includes nonlinearity and three dimensionality as well as 
rotational! ty. In contrast to this is wave induced particle motion 
which, while being random and three dimensional, is approximately ir- 
rotational and nondissipative (Lumley and Panofsky [5]). Figure 13 
clearly establishes the irrotationality of the wave induced velocities 
by showing an almost constant phase difference of 90^ between u and w 
components in the region of significant energy density. Tlic 90^ phase 
difference is predicted by solving the linear first order wave equation 
in terms of velocity potential. A qualifying condition for a velocity 
potential is the existence of irrotationality in the fluid. 

The horizontal component of water particle velocity, u, will not be 
considered in terms of absolute or quantitative values as the velocity 
meter measured only one component of the horizontal velocity during 



The measured energy-density spectrum of the vertical water particle 
velocity has been plotted in Figure 14. Figure 14 also shows the cal- 
culated values of the vertical velocity energy-density spectrum predicted 
using linear wave theory (Bordy [ 2 ]): 



The term in brackets is the first order vertical velocity amplitude - 
response function squared. T| is a measure of wave height as a function 



Run Five. 



>-.(2 r f sinh k(hfz) 
w(f) = i.u ' 



' L 2rr sinh kh 
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* Phase and Coherence Spectra Represent Wave Height versus W-measured. 




FIGURE 1^. MEASURED AND THEORETICAL VERTICAL VELOCITIES 
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of frequency. The measured velocity was equal to the theoretical value 
for the region of significant energy-density . Close examination of tlie 
spectra indicates that the measured and theoretical vertical velocities 
begin to diverge at a frequency of about 0.2 Hz. This marks the begin- 
ning of the transition regime to turbulence. At 0.3 Hz turbulence is the 
dominant process in the energy density spectrum. 

Spectra of other parameters compared to waves are shown in Figures 
10, 11, 12, and 15. The contribution to the energy-density spectra of 
each by wave induced particle motion in the frequency range .05 Hz to 
.22 Hz approaches 957, of the total energy witli the exception of tempera- 
ture . 

The phase and coherence spectra of the wave height versus modula- 
tion, shows the coherence approaching unity. This strong linear rela- 
tionship of v;aves and sound modulation was attributed to bubbles in an 
earlier section. The values of coherence are less for salinity, sound 
velocity, and temperature, with the least coherence for wave height ver- 
sus temperature. This could indicate that the temperature scalar field 
is more affected by turbulence and theirmal dif f usivi ty*, or it could in- 
dicate an instrumental inadequacy. Regardless, it is noted that the 
temperature has a significantly less linear relationship with wave height 
than the other parameters and accordingly only about 407o of its energy 
comes in the wave dominant portion of the spectrum. 

* The thermal diffusivity is approximately 10 times greater than salt 
diffusivity in the ocean. Thus temperature gradients are expected to 
be erased much quicker than salinity gradients. 
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VI. CONCLUSIONS 



A. The Fourier high pass filter proved to be an effective method of 
eliminating low frequency components from the signal record without al- 
tering phase, coherence, or energy relationships in the region of interest 

B. The spatial correlation distances for temperature determined for 
measurements at a point were established to be about 25 centimeters, 

C. Future experiments attempting to verify Wilson’s equation by summing 
the variances of its constituents (T,s,P) should require the instruments 
to be placed within a correlation distance; this insures that the local 
mechanisms affecting temperature and salinity inhomogenci ties arc iden- 
tical and will give a valid covariance term. Further, better salinity 
instrumentation is required for the time scales considered in the near 
surface region. 

D. Sound modulation of a 60 kHz sound beam traversing a 2 meter path 
length underwent a mean attenuation of about 57o. The correlation dis- 
tance of an effective inhomogeneity due to the combined temperature, salin 
ity, and bubble effects was found to be 8.8 centimeters by applying 
Mintzer and Stone’s Equation. 

E. The change in the index of the speed of sound at 60 kHz (a resonant 
bubble region) was found to be 997o induced by bubbles. 

F. The energy density spectra of all parameters were dominated by wave 
induced motion at 7.3 meters depth. 
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G. High coherence is found betv;een v/ave induced water particle motion 
and sound modulation and decreasing coherence between water particle 
motion and salinity, sound velocity, and temperature. 

H. First order wave theory accurately predicts orbital water particle 
velocities . 

I. Wave induced water particle velocities are irrotat ional . 
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